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ABSTRACT 

We present the first calculation of the kinetic Sunyaev-Zel'dovich (kSZ) effect due to the inhomo- 
geneous reionization of the universe based on detailed large-scale radiative transfer simulations of 
rcionization. The resulting sky power spectra peak at £ = 2000 — 8000 with maximum values of 
[£(£ + X)Ci /(27r)] max ~ 1 x 10~ 12 . The peak scale is determined by the typical size of the ionized 
regions and roughly corresponds to the ionized bubble sizes observed in our simulations, ^5 — 20 Mpc. 
The kSZ anisotropy signal from reionization dominates the primary CMB signal above £ = 3000. This 
predicted kSZ signal at arcminute scales is sufficiently strong to be detectable by upcoming experi- 
ments, like the Atacama Cosmology Telescope and South Pole Telescope which are expected to have 
~ 1' resolution and ~ /iK sensitivity. The extended and patchy nature of the reionization process 
results in a boost of the peak signal in power by approximately one order of magnitude compared 
to a uniform reionization scenario, while roughly tripling the signal compared with that based upon 
the assumption of gradual but spatially uniform reionization. At large scales the patchy kSZ signal 
depends largely on the ionizing source efficiencies and the large-scale velocity fields: sources which 
produce photons more efficiently yield correspondingly higher signals. The introduction of sub-grid 
gas clumping in the radiative transfer simulations produces significantly more power at small scales, 
and more non-Gaussian features, but has little effect at large scales. The patchy nature of the reioniza- 
tion process roughly doubles the total observed kSZ signal for I ~ 3000 — 10 4 compared to non-patchy 
scenarios with the same total electron-scattering optical depth. 

Subject headings: radiative transfer — cosmology: theory — cosmic microwave background — intergalac- 
tic medium — large-scale structure of universe — radio lines 



1. INTRODUCTION 

The secondary anisotropics of the Cosmic Mi- 
crowave Background (CMB) are emerging as one of 
the most powerful tools in cosmology. The small-scale 
anisotropics are probes of cosmological structures and 
are thus a valuable tool for studying their formation 
and properties. The key effect generating small-scale 
CMB anisotropies is the Sun yaev-Zel'dovich (SZ) effect 
ijZeldovich fc Sunvaevl 11969ft . produced by Compton 
scattering of the CMB photons on moving free electrons. 
When this is due to thermal motions the effect is referred 
to as thermal SZ effect (tSZ), while when it is due to 
the electrons moving with a net b ulk peculiar velocity it 
is ca lled kinetic SZ effect (kSZ) ijSunvaev fc Zeldovicbl 
1980). The former, coming largely from hot gas in 
low-redshift galaxy clusters, is the dominant effect, 
but has a different spectrum than the CMB primary 
anisotropies. The tSZ spectrum has a characteristic 
zero at ~ 217 MHz, so its contribution can in principle 
be separated. By contrast, kSZ anisotropies have a 
spectrum identical with that of the primary anisotropies. 
Separation of its effects is possible only by virtue of 
its different spatial structure. We consider the kSZ 
fluctuations as composed of two basic contributions, one 
coming from inhomogeneous reionization and the other 
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from the fully-ionized gas after reionization. The calcu- 
lation of the patchy reionization contribution to the kSZ 
signal based on the first large-scale radiative transfer 
simulations of rcionization is the main focus of this pa- 
per. This signal provides a signature of the character of 
rcionization that will complement other approaches like 
observations of the redshifted 21-cm line of hydrogen and 
surveys of high-z Ly-q emitters (e.g.lScott fe Reeslll990l 
i Tozzi et all2000Ulliey et al J20021 12003t iCiardi fc Madav. 
20031 IRhoads et alJ I2003t IStanwav fc et all 12004 
Gnedin fe Pradal 120041: iFurlanetto et alJ 12004a 
Mellema et al l l2006bt iMalhotra fe Rhoadsl 12006 
Shapiro et alJl2006t IBunker et al.H2006UFurlanetto et alJ 
120061) . In addition to better understanding of the 
physics of reionization, such studies could provide better 
constraints on the fundamental cosmological parame- 
ters, in particular on the primordial power spectrum of 
density fluctuations at smaller scales than are currently 
available. 

The kSZ from a fully-ionized medium has been stud- 
ied previously by both analytical and numerical means. 
When linear perturbation theory is used, the effect 
is associated with the quadratic nonlinearities in the 
electron density current an d is usually referred to as 
the Ostriker-V ishniac effect ijOstriker fe Vishniadll986t 
iVishniacI 11987). Calculating the full, nonlinear effect is 



much more difficult to do anal ytically, a lthough some 
models have b een proposed (e.g. IJaffe fc Kam ionkowskl 
Il998l 15^120001: IMa fc Frvll20"r)l IZhang et al.H2004fl . This 
effect involves coupling of very large-scale to quite small- 
scale density fluctuation modes and thus requires a very 
large dynamic range in order to be simulated correctly. 
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Current simulations are i n rough agreement, but have not 
quite converged yet (e.g . ISpringel et alJl2~00H IMa fc Frvl 
l2002tlZhang et aLll2004|) . All simulations predict signif- 
icant enhancement of the small-scale anisotropies com- 
pared to the analytical theory, due to the nonlinear evo- 
lution. 

Calculations of the kSZ contribution from patchy 
rcionization are even more challenging. In addition to 
the above difficulties, such calculations also require de- 
tailed modelling of the radiative transfer of ionizing ra- 
diation to derive the sizes and distributions of H II 
regions in space and time (the reionization geometry) 
and how these correlate with the velocity and density 
fields. To date this problem has been studied by only a 
few recent works. Most of these estimates were done 
by semi-analytical models, which used various simpli- 
fied approaches t o mode l the inhomog eneous reionization 
iGruzinov fc Hul 119981 iSantos et all I200& iZahn et all 
l2005H McQuin net al.H2005|) . While such models are use- 
ful since they are cheaper to calculate than full radiative 
transfer simulations, and so allow for fast exploration 
of the various scenarios, their results must be checked 
against full detailed simulations in order to ascertain 
their reliabilit y. The only two exist ing numerical studies 
of this effect. iGnedin fc Jaffc (200JJ) andlSalvatcrr a et al.l 
(2005) used fairly small computational boxes (4h _1 Mpc 
and 20h _1 Mpc, respectively). As a consequence, their 
results significantly underestimate the kSZ signal, as we 
discuss below. 

There are several upcoming experiments which will 
search for kSZ effect, in particular the Atacama Cosmol- 
ogy Telescope (ACT) 4 and South Pole Telescope (SPT) 5 , 
both of which are expected to be operational by early 
2007. These experiments will have ~ 1' resolution and 
~ fiK sensitivities, which should be sufficient to detect 
the kSZ signal. 

In this paper we present the first calculations of the 
kSZ effect from patchy reionization based on large-scale 
radiative transfer simulations. We simulate reionization 
using a (100 h _1 Mpc) 3 simulation volume, which is suf- 
ficient to capture the relevant large-scale density and 
velocity perturbations, an important improvement over 
previous efforts We use the ray-tracing code C 2 -Ray 
ijMellema et al.ll2006al) to follow the radiation from all 
ionizing sources in that volume identified with the re- 
solved halos, which are of dwarf galaxy size or larger. 
The halos and underlying density field are provided by 
a very large N- bo dy simulation with the code PMFAST 
(M erz et al.ll2005|) . The results of these simulations on 
the reionization character, geometry and observability of 
the reds hifted 21-cm line o f hydrogen we re discussed in 
detail i n llliev et aD l)2006bl Paper I) and iMellema et all 
(2006bl Paper II). A pr eliminary version o f our current 
results was presented in llliev et al.l l)2006dlh 

2. SIMULATIONS 

Our simulations follow the evolution of a comoving 
simulation volume of (100 /i _1 Mpc) 3 , corresponding to 
an angular size ~ 1 deg on the sky at the relevant red- 
shift range. Our methodology and simulation parameters 
were described in detail in Papers I and II. Here we pro- 



vide just a brief summary. We start with performing a 
very large pure dark matter simulation of early structure 
formation, with 1624 3 w 4.3 billion p articles and 3248 3 
grid cells 6 using the code PMFAST ijMerz et alJl2005|) . 
This allows us to reliably identify (with 100 particles or 
more per halo) all halos with masses 2.5 x 10 9 M© or 
larger. We find and save the halo catalogues, which con- 
tain the halo positions, masses and detailed properties, in 
up to 100 time slices starting from high redshift (z ~ 30) 
to the observed end of reionization at z ~ 6. We also 
save the corresponding density and bulk peculiar veloc- 
ity fields at the resolution of the radiative transfer grid. 
Unfortunately, radiative transfer simulations at the full 
grid size of our N-body computations are impractical on 
current computer hardware, thus we solve for the radia- 
tive transfer on coarser grids, of sizes 203 3 = (3248/16) 3 
or 406 3 . These grid resolutions allow us to derive reliably 
the angular sky power spectra for I as 430 — 90,000 for 
203 3 , and up to I ~ 180,000 for 406 3 . Throughout this 
study we assume a flat ACDM cosmology with parame- 
ters (i] m ,fi A ,ft,/i (T 8 ,)i) = (0.27,0.73,0.044,0.7,0.9,1) 
l)Spergel et alJ l2003k where fl m , Qa, and fib are the to- 
tal matter, vacuum, and baryonic densities in units of 
the critical density, h is the Hubble constant in units of 
100 kms _1 Mpc _1 , erg is the standard deviation of linear 
density fluctuations at present on the scale of 8h _1 Mpc, 
and n is the index of the primordial power spectrum. We 
use th e CMBfast transfer function ijSeliak fc Zaldarriagal 
1996). All calculations are done in flat-sky approxima- 
tion, which is appropriate for our relatively small angular 
sizes. 

All identified halos are assumed to be sources of ion- 
izing radiation and each is assigned a photon emissivity 
proportional to its total mass, M, according to 
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(1) 



where t s is the source lifetime, m p is the proton mass, \i 
is the mean molecular weight and / 7 is a photon produc- 
tion efficiency which includes the number of photons pro- 
duced per stellar atom, the star formation efficiency (i.e. 
what fraction of the baryons are converted into stars) 
and the escape fraction (i.e. how many of the produced 
ionizing photons escape the halos and are available to 
ionize the IGM). 

The radiative transfer is followed using our fast and ac- 
curate ray-tracing photoionization and non-equilibrium 
chemistry code C 2 -Ray. The code has been tested in 
detail for correctness and accuracy against available an- 
alytical solutions and a number of other cosmological ra- 
diative transfer codes ((Mellema et alJl2006at llliev et all 
l2006a|L The radiation is traced from every source on the 
grid to every cell. 

We have performed four radiative transfer simulations. 
These share the source lists and density fields given by 
the underlying N-body simulation, but adopt different 
assumptions about the source efficiencies and the sub- 
grid density fluctuations. The runs and notation are 
the same as in Paper II: runs f2000 and f250 assume 

6 3248 = N„ odea X (512 - 2 X 24), where N nodes = 7 (with 4 
processors each). 512 cells is the Fourier transform size and 24 cells 

Ihttp://www.hep.upenn.eduAmskip-\thiimuskip\sim\msk^ force matching on each side of 

http://spt .uchicago.edu/science/index.html the cube. 
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/ 7 = 2000 and 250, respectively, and no sub-grid gas 
clumping, while f2000C and f250C adopt the same re- 
spective efficiencies, / 7 = 2000 and 250, but also add a 
sub-grid gas clumping, C{z) = (n 2 )/(n) 2 , which evolves 
with redshift according to 

C subgrid (z) = 27.466e-°- 114z+0 - 001328 * 2 . (2) 

The last fit was obtained from another high- 
resolution PMFAST N-body simulation, with box size 
(3.5 h _1 Mpc) 3 and a computational mesh and number 
of particles of 3248 3 and 1624 3 , respectively. These pa- 
rameters correspond to a particle mass of 10 3 M Q and 
minimum resolved halo mass of 1O 5 M0. This box size 
was chosen so as to resolve the scales most relevant to 
the gas clumping - on scales smaller than these the gas 
fluctuations would be below the Jeans scale, while on 
larger scales the density fluctuations arc already present 
in our computational density fields and should not be 
doubly-counted. The expression in equation J3J excludes 
the matter inside collapsed minihalos (halos which are 
too small to cool atomically, and thus have inefficient 
star formation) since these are shielded, unlike the gen- 
erally optically-thin 1GM. This self-shielding results in a 
lower contribution of the minihalos to the total number 
of recombinations than one would infer from a simple 
gas clum ping argument ijShapiro et al.ll2004l llliev et aD 
2005b -J). The effect of mini halos could be inclu ded as 
sub-grid physics as well, see iCiardi et alJ ((2006). This 
results in slower propagation of the ionization fronts and 
further delay of the final overlap. The halos that can cool 
atomically are assumed here to be ionizing sources and 
their recombinations are thus implicitly included in the 
photon production efficiency / 7 through the correspond- 
ing escape fraction. 

3. KSZ FROM PATCHY REIONIZATION 

The kSZ effect is the CMB temperature anisotropy 
along a line-of-sight (LOS) defined by a unit vector n 
induced by Thomson scattering from flowing electrons: 

AT f 

— = dr)e~ Tee Wan e a T n-v, (3) 

J- CMB J 

where r/ — J Q dt' /a(t') is the conformal time, a is the scale 
factor, ctt = 6.65 x 10 -25 cm -2 is the Thomson scattering 
cross-section, and r es is the corresponding optical depth. 

We calculate the kSZ anisotropy signal from our simu- 
lation data as follows. We first calculate the line-of-sight 
integral in equation for each individual output time- 
slice of the radiative transfer simulation and for all LOS 
along each of the three box axes. The contribution to the 
total LOS integral from each light-crossing time of the 
box is then obtained by linear interpolation between the 
nearest results from the simulation output times. The 
individual light-crossing time contributions are then all 
added together to obtain the full LOS integral given in 
equation All these integrals are done for each LOS 
through the box along the direction of light propagation, 
which allows us to produce kSZ maps, in addition to the 
statistical signals. In order to avoid artificial amplifica- 
tion of the fluctuations resulting from repeating the same 
structures along the line of sight, after each light-crossing 
time we randomly shift the box in the directions perpen- 
dicular to the LOS and rotate the box, so that the LOS 



cycles the directions along the x, y and z axes of the 
simulation volume. 

For comparison we also consider two simplified cases 
in addition to our simulations, the cases of instant and 
of uniform reionization. We define instant reionization 
as a sharp transition from completely neutral to fully- 
ionized IGM at redshift 2; insta nt- We pick z insta nt = 13, 
which yields the same integrated electron scattering op- 
tical depth as our simulation f250. We also define a "uni- 
form reionization" scenario to be one that has the same 
time-dependent reionization history (and hence, also the 
same r es ) as simulation f250, but spatially uniform (i.e. 
not patchy). We then derive the kSZ temperature fluctu- 
ations for these two scenarios using the same density and 
velocity data and same procedures as for the actual sim- 
ulations. We consider these simplified models in order 
to demonstrate the effects of reionization being extended 
and patchy in nature. 

4. RESULTS 
4.1. kSZ maps 

We show the kSZ maps of temperature fluctuations, 
ST/Tcmb, yielded by all of our cases in Figure ^ These 
have a total angular size of approximately 50' x 50', cor- 
responding to our computational box size. The resolu- 
tion of the maps is that of the full radiative transfer grid 
(203 x 203 pixels), corresponding to pixel resolution of 
~ 0.25'. All maps utilize the same color map in order 
to facilitate their direct comparison. The maps deriving 
from our simulations of inhomogeneous reionization (Fig- 
ure ^ top and middle) all show fairly strong fluctuations, 
both positive and negative, of order ST ~ 10 fiK at angu- 
lar scales of a few arcminutes (~ 10 h _1 Mpc) and up to 
ST ~ 20 iiK at the full map resolution. The fluctuations 
are at somewhat smaller scales when the sources are less 
efficient photon producers (f250 and f250C), compared to 
the high-efficiency cases (f2000 and f2000C). The varia- 
tions are noticeably enhanced when sub-grid gas clump- 
ing is included in the reionization model (f2000C and 
f250C) and there are a number of regions with very 
strong features. In comparison, the artificial models of 
instant and uniform reionization (Figure^ bottom) show 
fluctuations with much lower amplitudes and with a less 
well-defined typical scale. Since these simple scenarios 
were constructed to produce the same total electron scat- 
tering optical depth as simulation f250 and they use the 
same density and velocity fields as the simulations, any 
observed differences are due to the patchiness of reioniza- 
tion. The lower signal in these last two cases is expected 
since the kSZ effect in uniformly-ionized g as exactly can - 
cels in the first order of the linear theory l)Kaiseil fl984'l. 
This cancellation is broken when patchiness is present, 
however, which enhances the anisotropies. 

The mean and rms of the temperature fluctuations for 
all cases are summarized in Tabic ^ Based on these, 
we see that, indeed, the uniform and instant reionization 
cases have means very close to zero, as expected, while 
the realistic, patchy reionization simulations yield mean 
temperature fluctuations that are low, of order ~ 10~ 7 , 
but not zero. The rms values are ~ 10 -6 for all patchy 
reionization cases and lower than that by factor of ~ 2 
(~ 3) for the instant (uniform) reionization cases. 

These observations based on the maps are confirmed 
by their corresponding pixel PDF distributions, shown 
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Fig. 1. — kSZ maps from simulations: f2000 (top left), f250 (top right), f2000C (middle left), f250C (middle right) instant reionization 
at ^instant = 13 which gives the same total electron scattering optical depth as simulation f250 (bottom left), and spatially- uniform 
reionization with the same reionization history and thus same total electron scattering optical depth as simulation f250 (bottom right). 
(Images produced using the Ifrit visualization package of N. Gnedin). 



in Figures |21 and |3 On each panel we also plot the Gaus- 
sian distribution with the same mean and standard de- 
viation. All distributions are surprisingly close to Gaus- 
sian around their mean values given the maps. However 
noticeable departures from Gaussianity do occur in the 



wings of the PDFs. In particular, the PDFs for patchy 
reionization with sub-grid clumping (f2000C, and f250C) 
are significantly non-Gaussian. For the realization being 
considered here there is an over-abundance of bright re- 
gions by up to an order of magnitude compared to the 
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TABLE 1 

Mean and rms values for <5T k sz/7cMB- 



simulation f2000 f2000C f250 f250C uniform instant 

mean 8.86 x 10~ 8 9.88 x 1CT 8 -2.17 x 10~ 7 5.57 x 10" 8 -3.29 x 10" 12 -3.37 x 10" 
1.13 x 1CT 6 1.22 x 10~ 6 1.06 x 10~ 6 1.09 x 10~ 6 4.24 x 10~ 7 6.67 x 10~ 7 
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Fig. 2. — PDF distribution of <5Tksz/TcMB (solid) vs. Gaussian 
distribution with the same mean and width (dotted) for our four 
simulations, as labelled. 



corresponding Gaussian distributions. The reionization 
scenarios without sub-grid clumping show much weaker 
non-Gaussian features there. This indicates that the ob- 
served PDF of kSZ from patchy reionization may give 
us important information on the level of small-scale gas 
dumpiness during reionization. The PDFs derived from 
the uniform and instant reionization scenarios are sig- 
nificantly less wide than the simulated ones (which was 
also shown by their lower rms values, as was discussed 
above), and are also much closer to Gaussian. 

4.2. Sky power spectra 

In Figure 0] we show the 2D sky power spectra de- 
rived from our simulation data. These were obtained 
by calculating the sky power spectra for each box light- 
crossing time and adding these, which yields a smoother 
final result than would calculating the power spectrum 
directly from the final map. The kSZ anisotropy sig- 
nal from inhomogeneous reionization dominates the pri- 
mary CMB anisotropy above I = 3000. The sky power 
spectra peak strongly at t = 3000 — 5000, to a maxi- 
mum of [£(£ + l)CV/(27r)] max - 7 x 10" 13 when the ioniz- 
ing sources are highly- efficient (cases f2000 and f2000C). 
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Fig. 3. — PDF distribution of (5T k s z /Tc M B (solid) vs. Gaussian 
distribution with the same mean and width (dotted) for the fiducial 
scenarios of uniform (left) and instant reionization (right). 



The peak values for the simulations with lower-efficiency 
sources (f250 and f250C) are only slightly lower, at [£(£ + 
l)C£/(27r)] max ~4x 10~ 13 , and the peaks are somewhat 
broader and moved to smaller scales, I ~ 3000 — 7000, 
than for the high-efficiency cases. The introduction of 
sub-grid gas clumping in the radiative transfer simula- 
tions produces significantly more power at small scales 
and slightly broader peaks when compared to the simi- 
lar cases without sub-grid clumping. The scales at which 
the power spectra peak roughly correspond to the typical 
ionized bubble sizes observed in our simulations, namely 
^5 — 20 Mpc comoving. These typical sizes depend 
on the assumed source efficiencies and sub-grid clump- 
ing. At large scales the patchy kSZ signal depends only 
on the source efficiencies, being higher for more efficient 
photon emitters, since these tend to produce larger ion- 
ized regions on average. 

In contrast, the uniform reionization scenario (with the 
mean reionization history of simulation f250) yields a 
kSZ signature which is much lower than the non-uniform 
reionization scenarios, indicating a very large boost of 
the signal due to the effects of patchiness. The boost is 
largest, approximately one order of magnitude, at and 
above the typical scale of the patches (£ < 7000), but it 
still exists at smaller scales, where it is a factor of two 
or more. The other simplified scenario, of instant reion- 
ization with the same integrated optical depth as f250, 
produces larger kSZ anisotropy than a uniform reion- 
ization does. However, it is still well below the realistic 
patchy reionization signals, by factor of ~ 3 for £ < 7000. 
At smaller scales the reionization signal from the instant 
reionization scenario becomes similar to the ones from 
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Fig. 4. — Sky power spectra of <57kgz/TcMB fluctuations re- 
sulting from our simulations: f2000 (black), f250 (blue), f2000C 
(red) and f250C (green). For comparison, we also show the results 
from simple models which utilize the same density and velocity 
fields as the actual simulations, but different assumptions about 
the gas ionization in space and time, a uniform reionization with 
the same reionization history as simulation f250 (dashed, dark red) 
and an instant reionization model with the same integrated optical 
depth r es as simulation f250 (dashed, cyan). The primary CMB 
anisotropy signal is also shown (dotted, brown). 



simulated non-uniform reionization. The distribution of 
the power in these two simplified scenarios is much flat- 
ter, with less indication of a characteristic scale. This 
implies that the sharp peaks yielded by the inhomoge- 
neous reionization scenarios are dictated by the size of 
the ionized patches, since both the density and the ve- 
locity fields are shared among our models. This point is 
discussed further in § 14.31 

This behavior can be understood further by consider- 
ing the contributions from different redshift intervals to 
the integrated signal. In Figure we show the contribu- 
tions to the total signal from different redshift intervals 
for simulation f250, as well as the uniform and instant 
reionization scenarios. We plot the contribution from 
every three light-crossing times of the box, correspond- 
ing roughly to z > 15, 15 > z > 11 and z < 11. The 
mass- weighted ionized fraction in the highest redshift in- 
terval is x rn < 0.01, and consequently its contribution to 
the kSZ effect (case f250, bottom curve) is low, with a 
maximum of ~ 10 -14 . This high-redshift contribution is 
strongly peaked at very small scales (£ > 10 4 ), reflecting 
the fact that at that time the typical ionized bubble is 
still fairly small, less than a few Mpc in size. The con- 
tribution from the middle redshift interval, 11 < z < 15, 
(corresponding to 0.8 > x m > 0.01) peaks at roughly the 
same scales as the total signal. This interval contributes 
about half of the total signal around the peak, but a 
smaller fraction at scales above and below that. A half 
or more of the integrated kSZ signal at all scales is con- 



tributed by the lowest redshifts (z < 11, x m > 0.8). The 
low- redshift power spectrum is fairly flat, with a weak 
peak at relatively large scales, I ~ 2000. This late-time 
contribution resembles the results for the two homoge- 
neous scenarios, but has higher amplitude and a different 
shape. This indicates that even at late times, when most 
of the IGM is already ionized, the patchiness still plays 
an important role in assembling the kSZ signal. 

In the homogeneous (uniform and instant) reionization 
cases the shape of the power spectra is essentially the 
same for all redshift intervals, quite flat, with a broad 
peak at t ~ 5000. The signal from the uniform reion- 
ization scenario is completely dominated by the lowest 
redshifts, with less than a few per cent coming from the 
middle redshift interval, and essentially no contribution 
from high redshifts. Since the velocity and density fields 
are shared, the instant reionization scenario contribu- 
tions differ from the corresponding ones from the uniform 
scenario only by being weighted by the ionized fraction 
in the uniform case. The instant reionization scenario as- 
sumes full ionization after Zi nsta nt = 13 and fully-neutral 
gas before that. Thus, as a consequence to the kSZ sig- 
nal being stronger at later times, the instant reioniza- 
tion integrated signal is higher as well. For comparison, 
we also show the corresponding red shift-binned results 
from the model of iZhang et all i|200l (see against 
our results from the instant reionization scenario. The 
two results agree well at small scales, in the potentially- 
observable range. 

4.3. Linear theory and large-scale velocity field 
corrections 

In linear theory the velocity and density perturbations 
at redshift z are related through the continuity equation: 



6(M) = ^^u(M), 



(4) 



k 2 D(t) 

where a = (1 + z)^ 1 is the scale factor, and D(t) is the 
growth factor of linear perturbations, which for a flat 
ACDM background cosmology is given by 



D(z) 



5Q, E(z) 



1 



[E{z')\ 



-dz . 



(5) 



with E(z) = H(z)/H = [fi (l + zf + ^a] 1/2 , where 
H(z) and Hq are the values of the Hubble constant at 
redshift z and at present, respectively. In terms of the 
power spectra of the density, Pss, and the velocity, P VVl 
equation can be re-written as 

~D(t)" 2 



or 



A 2 v (k,t) 
where we defined 

and 



a 2 (t) 
k 2 



a 2 (t) 
k 2 



D{t) 

Pit) 
D(t) 



Pss{k,t), 



A?(A,i), 



A 5 = 
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Fig. 5. — Contributions to the total kSZ signal (blue) from different rcdshifts for case f250 (left), uniform (middle) and instant reionization 
(right). Shown are the signals for every three box light crossing times, roughly corresponding to (bottom to top on the left) redshifts z > 15 
(green), 15 > z > 11 (red), z < 11 (magenta). In the last two cases the z > 15 contribution is very low (uniform) or zero (instant) and 
thus n ot sho w. F or comparison, on the last panel we also show the corresponding redshift-binned results from the model of [Zhang et alj 
Pool) (see S lOt . 
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Fig. 6. — Linear-theory power spectra at redshift z = 13.6 
of the density field (dotted, black) and velocity field (in kms -1 ; 
solid, black) compared to the simulation results for the density 
(long-dashed, red) and velocity (short-dashed, red). 



which would be our notation for the power spectra here- 
after. The quantities A 2 and are the power per loga- 
rithmic interval of the wave-number k of the density and 
velocity fields, respectively. 
Using equation (J5J it is straightforward to show that 



D 
D 



a a 5f2o o, (1 + z) 2 
1 ~ a + ~a[E(z)] 2 D{z) 
3H n (l + zf 



E{z) 



1 



3(l + z)£>(z) 



(10) 



Substituting equation (fTU|) into equation Q we finally 
find 

A 2 v (k,z) = 



Ag(fc,z) 9H 2 n 2 (l + z) 
k 2 AE 2 {z) 



1 



3(1 + z)D(z) 



.(11) 



In Figure we show the CMBfast density power spec- 
trum at redshift z = 13.6, and the corresponding linear- 
theory velocity power spectrum given by equation (fill) , 
along with the density and velocity power spectra ob- 
tained from our N-body simulation data (re-gridded to 
our radiative transfer grid of 203 3 ). Comparison of 
the two sets of data shows a good agreement between 
them and indicates that at these scales and redshifts our 
density and velocity fields are clearly still in the linear 
regime. Note, however, that the bulk velocities derived 
from the simulation data do include the non-linear effects 
at smaller scales, up to the full resolution of the under- 
lying N-body simulation, at 3248 3 cells, or ~ 30/i _1 kpc 
comoving per cell, averaged at the radiative transfer grid 
resolution. The velocity power spectra were derived from 
the simulated data by first calculating the power spec- 
trum of each of the three velocity components and then 
averaging these, which minimizes the variance at large 
scales. 

The velocity power spectrum has a broad peak at fairly 
large scales, k ~ 0.01 — O.lhMpc -1 , corresponding to 
scales ~ 60 — 600 Mpc h _1 . Thus, our simulation volume 
is sufficiently large to reach this velocity peak, but is still 
missing some velocity power from the largest scales. We 
estimate the missing power as follows. The rms of the 
velocity field is given by 



v 2 

rms 



= J A 2 v (k)dln(k). 



(12) 



Thus, integrating over the full linear power spectrum 
yields the total power in the velocity field, v 2 ms tot in 
the linear theory. A finite simulation box of size Lbox 
would not include the modes with wave-numbers below 
fcmin = 27r/Lbox- Integrating A 2 over the wave-numbers 
k > fcmin, we obtain the total velocity power for a given 
box size. The results are shown in Figure where we 
plot the fraction of total linear-theory velocity power 
u r 2 ms box/^rms tot present in a simulation box vs. the box 
size I>box- A simulation of volume (100 h _1 Mpc) 3 would 
thus be missing about ~ 50% of the total v 2 ms velocity 
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Fig. 7. — Fraction of total linear-theory velocity power 
u ?ms box/^ms.tot present in a simulation box of size L box vs. L box . 



power as given by the linear theory. Simulations with 
smaller volumes would miss much larger fraction of the 
velocity power, ~ 70% for 50h _1 Mpc box, ~ 90% for 
20tr 1 Mpc box, and over 99% for 4h~ 1 Mpc box. 

The large-scale bulk motions not included in our simu- 
lation occur on scales of ~ 100 Mpc or larger, well above 
the cha racteristic size of the ionized patches (~ 10 Mpc 
or less) l)Iliev et al.l2006bl see also FigureEJ, thus on such 
scales the reionization patchiness averages out and the 
ionization fraction approaches the mean for the universe. 
We can approximately account for the missing large-scale 
velocity power as follows. At each light-crossing time we 
can assume that the whole simulation volume is moving 
with some random velocity, Vhox- In order to account for 
the missing velocity power, the value of the component of 
this random velocity at redshift z along the line-of-sight 
is given by 



Wbox(z) = V. 



rms, missing 



where v, 



rms, missing 



(-21n<7) 1/2 cos(27r(9), (13) 

[Wrms,tot(^) 2 - Vrms,box(.2)] 1/2 , and 

q and 9 are uniformly-distributed random variables be- 
tween and 1. The form of equation (|13|l ensures that 
VboA z ) lS a Gaussian-distributed r andom variable with a 
zero mean and rms of u r ms,missing l|Box &: Mullerl [l958). 
In terms of the contribution from this redshift to the 
temperature anisotropy this yields 



V TqMB ) , 



AT 
V Tqmb J 



(z) + T cs (z) 



Vbo X ( Z ) 



1 ) OX 



(14) 

where t cs (z) is the corresponding contribution to the 
total electron-scattering optical depth. Note that for 
uniform reionization the effects from large-scale veloc- 
ity fields wou ld exactly cancel, since the flow is potential 
ljKaiserlll984|) . However, the patchiness breaks this can- 
cellation and the large-scale velocities increase the signal 



Fig. 8. — Sky power spectra of <57kgz/TcMB fluctuations from 
the epoch of reionization based on our simulations, f2000 (black), 
f250 (blue), f2000C (red) and f250C (green), corrected for the miss- 
ing velocity field power compared to the after-reionization kSZ sig- 
nals (assuming overlap at z ov = 8): linear Ostrikcr-Vishniac ef- 
fect, labelled 'OV (long-dashed, pink), the same, but using the 
nonlinear power spectrum of density fluctuations, labelled 'NL' 
(long-dashed, dark red) and a fully-nonlinea r model mat c hed t o 
high-resolution hydrodynamic simulations of Zhang ct al. (2004), 
labelled 'Zhang et al' (short-dashed, dark green). The primary 
CMB anisotropy signal is also shown (brown, dotted). 



by increasing the ionized bubble velocities along the LOS. 
The ionized regions are at much smaller scale than these 
large-scale velocities and uncorrelated with them, thus 
avoiding the usual cancellation. 

The power spectra with large-scale velocity correc- 
tions included are shown in Figure |S| For compari- 
son, we also show three representative calculations of the 
post-reionization contribution from fully-ionized gas af- 
ter rei onization: the quadratic-order Ostriker-Vishniac 
effect ljVishniadll987l) expressed in terms of a product 
of the linear power spectra; the same expression, but 
with the nonlinear density power spectrum substituted 
for one of the linear ones, which partially accounts for 
the nonli near effects; and th e recent detailed nonlinear 
model of Zha ng et alJ (J2004) . The last three models are 
rescaled to our ado pted cosmology using t{l + l)Cg cx erf 
( Zha ng et alJ l2004). and assume reionization overlap at 
z v = 8, to allow for direct comparison. Confronting first 
the three post-reionization models, we see that they agree 
fairly well at large, linear scales, but strongly diverge, 
by up to one order of magnitude at small scales, where 
the nonlinearities become very important. Including the 
corrections due to the nonlinear density power spectrum 
yields modestly higher kSZ signal than the linear OV cal- 
culation, and similar power spectrum shape overall. Both 
pea k at roughly the sa me scale, I ~ 2000. In contrast, 
the iZhang et al.1 l|2004D result finds much larger signal, 
especially at small scales, which peaks at I ~ 20, 000. 
Compared to our patchy reionization kSZ predictions, all 
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uniform-ionization power spectra are much less sharply 
peaked, i.e. they lack a well-established characteristic 
scale. In terms of their peak values, the OV and OV 
with nonlinear corrections models yield valu es lower than 
any of our patchy reionization results, while Iztiang et aD 
(2004) find a similar signal around the patchy power 
spectra peaks. At small scales (I > I0 4 ) the patchy 
reionization kSZ signals are similar to or lower than the 
OV result and are much lower than the expected full 
nonlinear post-reionization contribution. Thus, the best 
range to aim for when attempting to detect the patchy 
reionization signal is I = 3000 — 10 4 . 

4.4. Characteristic scales and LOS spectra 

The characteristic scale corresponding to the peak of 
the power spectra in Figure El is largely dictated by the 
typical scales of the ionized bubbles, since, as we dis- 
cussed in § 14.31 the velocity and density are in the linear 
regime and do not have characteristic scales within the 
range corresponding to our simulation volume. However, 
the ionized patches do have a typical scale. In Figure 
we show A^, the power spectrum of the ionized fraction 
(which is the same as the one of the neutral fraction) 
for all of our simulation cases at the redshifts when the 
volume ionized fractions are x v = 0.3, 0.5 and 0.8. In 



all cases the power spectra peak around wave-numbers 
k ~ 1 hMpc -1 , which scale corresponds to I ~ 7000. 
As reionization progresses, the peak moves gradually to 
slightly larger scales, resulting in a moderate widening of 
the peak of the kSZ temperature anisotropy sky power 
spectra and slightly less well-defined typical scale than 
seen in the single-redshift 3D power spectra. In order 
to quantify the power spectra evolution better, in Fig- 
ure we show the evolution vs. redshift of the peak 
amplitude, A x , max , the scale at which it occurs, fc max , 
and the effective width of the peak, given by 



width = k sinh 



1/2 



where 



5{\nk) 2 = / [(lnfc) 2 - (lnfc) 2 ]A^lnfc. 



Here 



k = exp(ln/c) = exp 



(\nk)A 2 x dlnk 



and 



A 2 x dln(k), 



(15) 

(16) 

(17) 
(18) 
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Fig. 10. — Evolution of A x . Plotted are the peak value, A x ,max, and the scale at which it occurs, fc ma x, its effective width, and the 
corresponding rms, x TulB for all simulations, as labelled. We also show the bias, b(z), at large scales of the ionized density with respect to 
the underlying density field (values indicated on the right axis) and k (see text for details). The values of fco/2, where ko is the value of 
the characteristic ionized bubble scale which give a good fit to the W(k) window function at large scales are also indicated for each case 
(black squares; see text for details). 



which gives the total integrated power per logarithmic 
interval. We also show k (in /i _1 Mpc) and x^ ms . The 
amplitude of the peak initially rises, peaking at the 
point of maximum patchiness, x rn ~ 0.5, and slightly 
decreases thereafter. The scale of the peak, fc max , (in 
/i _1 Mpc) steadily increases in time, as the typical ion- 
ized regions grow. Approximate linear fits (in the form 
a + bz) to k max (z) are given by a = —0.528 (-0.329,- 
1.135,-0.847) and b = 0.085 (0.089,0.17,0.19) for simula- 
tion f2000 (f2000C, f250, f250C). The effective width of 



the power spectra and k are largely constant and roughly 
equal to each other, except for simulation f250C at early 
times. Finally, the rms initially strongly rises, peaking 
at slightly later times than the patchiness itself peaks (at 
x m ks 0.6). 

During the later stages of reionization [x m > 0.3), 
which contribute most of the kSZ signal the ionized 
density fluctuations at large scales are proportional to 
the densit y fluct uations with a bias coefficient b(z) 
dJ2 006b. see also FigureEl Note however that 
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this is not the case at earlier times, when the ionized 
regions are relatively small and isolated and do not fol- 
low the large-scale density perturbations.). Based on this 
fact we can derive a handy approximate expression for 
the power spectrum of the ionized gas density, As.hii, 
as follows. Let 

&5,Hii(k, z) = b(z)W(k, z)A s (k, z), (19) 
where W(k, z) is some window function asymptoting to 1 
at large scales. We recognize that because of the zero or 
one nature of reionization (i.e. a region is either ionized 
or neutral when ionization fronts are unresolved, as is the 
case here), the random field 8n e /n e is more complex than 
this relation between the power spectra indicates. How- 
ever, it is reasonable expect that for low k, the electron 
density structure will have such a proportionality, with 
the window function being related to a quadratic super- 
position of form factors for the various H II patches. If 
the reionization contribution from dense regions is less- 
ened because of negative feedback from earlier ionizers 
on nearby galaxies, a situation we have ignored, the re- 
lation could be considerably modified. In spite of such 
complications, a simple relationship characterized by a 
characteristic scale related to the average size of the H II 
regions turns out to prevail as we now show. 

The bias b{z) is readily calculated from the simulation 
data (at the largest scale available from our simulation 
volume) and is shown in Figure ITUI (values shown on the 
right axis). The bias starts close to zero since very lit- 
tle of the gas is ionized, and has a peak value of 2-3, 
reached again around the time of maximum patchiness. 
For redshifts around the time of maximum patchiness 
(for mass-weighted ionized fraction 0.3 < x m < 0.8) and 
at large scales (k ;> 1 — 2) the window function W(k, z) 
is reasonably well-fit (within factor of ~ 2, and generally 
much better) by the following functional form: 

^ Z » = IT^ (20) 

where k${z) is a characteristic scale parameter which 
gives the best fit at the turnover. Samples of this window 



function, b(z)W(k,z), for simulation f250 are shown in 
Figure ITT1 At small scales this function turns over and 
the above fit is not applicable. The values of ko(z)/2 at 
x m = 0.3, 0.5 and 0.8 (in h~ 1 Mpc) are shown in FigurelTUl 
(squares). They are very close to the scale at which the 
ionized fraction power spectra peak, fc max . 

The addition of the missing large-scale velocity power 
outlined in the previous section boosts the sky power 
spectra at all scales, but especially at the largest ones, 
which further widens the peak. Nonetheless, the char- 
acteristic ionized bubble scale remains imprinted on the 
final results through their peak at £ ~ 2000 - 9000. 

To gain some further insight in the complex nature of 
the derived kSZ signal, in Figures El and El we show 
some sample line-of-sight (LOS) cuts through the sim- 
ulation volume (simulation f250) of the density field in 
units of the mean, nn/nii, ionized fraction, xu, elec- 
tron number density in units of the mean gas density, 
n e /nH = xnnn /fin, velocity component along that 
LOS, v x , and the accumulation of the kSZ tempera- 
ture anisotropy integral (from left to right), ATksz(< 
Early in the evolution (z = 13.62; Figure El the H II 
regions tend to be associated with high-density peaks, 
where the first halos (ionizing sources) form (inside-out 
reionization scenario). The ionized fraction distribution 
acts as a filter, picking the density field portions that 
would contribute to the CMB temperature anisotropies. 
The ionized regions are correlated with the high-density 
peaks (inside-out reionization), but the correlation is 
not perfect, there are density peaks which still remain 
neutral, while at the same time there are ionized re- 
gions which have a low overdensity, or are even under- 
dense. The velocity field fluctuations are dominated by 
the largest scales and are generally not closely correlated 
with the density field or the ionized fraction. The ve- 
locity field determines the sign of the contribution of a 
given ionized region to the total temperature anisotropy. 
A number of different situations can be observed in our 
data. For example, in Figure El (left) the ionized re- 
gions are relatively closely spaced together, separated by 
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Fig. 13. — LOS cuts through the simulation volume (simulation 
f250) at z = 11.310 (x v = 0.66, x m = 0.70). 



~ 10 Mpc, falling within the same large-scale velocity 
fluctuation. All H II regions except one are moving in 
the positive direction, resulting in a relatively large inte- 
grated temperature anisotropy of ~ 4.4 [iK. The other 
sample case from the same redshift (Figure l"H21 right) in- 
cludes three ionized regions along the LOS. These are 
separated by 40-50 comoving Mpc, and thus correspond 
to completely uncorrelated density fluctuations. The ve- 



locity field does not appear to correlate them, either, but 
nonetheless the temperature anisotropies they produce 
happen to almost exactly cancel each other, resulting in 
AT kSZ ~ along that LOS. At later times (Figure H3J) 
most of the volume is already ionized and each LOS en- 
counters a number of ionized regions, but ultimately al- 
most all of the total signal is contributed by the two 
highest-density peaks along the LOS, both of which here 
happen to move in the negative direction, while the con- 
tributions of the multiple other H II regions largely cancel 
each other. Once again, the density and velocity fields 
are not correlated. The ionized fraction is less correlated 
with the high-density regions than at early times. These 
examples demonstrate some of the complexities one has 
to deal with when trying to derive the kSZ signal from 
patchy reionization analytically, underscoring the need 
for large-scale simulations of this effect. 

5. COMPARISONS TO PREVIOUS WORK 

Previous si mulations of the kSZ e ff ect from patchy 
reionization (|Gnedin & Jaffd [200lt ISalvaterra et~all 
2005) predicted significantly lower kSZ signals than the 
ones we find. Their power spectra reach maximum values 
of [£(£+l)Ci/(2Tr)] max ~ 2xl0~ 14 and ~ 1.6 x 10~ 13 , re- 
spectively, compared to [£(£ + l)Ci/(2Tt)] max ~lx 10~ 12 
for our simulations. This discrepancy is due to the small 
volumes used in these simulations. As we showed in 
S 14.31 considering small volumes significantly reduces the 
power in the velocity field fluctuations. At the scale of 
the simulation box the bulk velocities are zero by def- 
inition, and the density is the mean one for the uni- 
verse, thus any larger-scale fluctuations are not included, 
and the ones close to the box size are underestimated. 
The reionization patchiness imprints its characteristic 
scale on the density and velocity fluctuations, effectively 
smoothing the small-scale fluctuations below the typical 
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bubble sizes. On the other hand, the large-scale fluc- 
tuations still should contribute to the kSZ signal, since 
the H II regions are moving with the large scale bulk 
motions. 

In addition to the missing large-scale power, these 
earlier reionization simulations also did not follow suf- 
ficient volumes to properly sample the size distribution 
of the ionized bubbles. The typical sizes of these ion- 
ized patches are of order 5-20 comoving Mpc, and be- 
come even larger at later times. These large character- 
istic bubble scales are result of the strong clustering of 
ionizing sources at high redshifts. Capturing the full bub- 
ble size distribu tion requires simul ation volumes of order 
(lOOh^Mpc) 3 (|Iliev et alJl2006b|) . and as a result the 
kSZ power spectra found by the smaller-box reionization 
simulations were largely flat, with no clear characteristic 
scale. 

Several semi-analytical models for calculating the kSZ 
signature of patchy reionization have been proposed in 
recent years l|Gruzinoy fc HiJ Il 998t ISantos et alj 120031: 
iMcQuinn et al.ll2005|) . iGruzinov fc Hul l)1998j) proposed 
a very simple model, whereby the ionized patches are 
randomly distributed and have a given characteristic size 
R. Furthermore, they assumed that the density, velocity 
and ionization fraction fields are all uncorrelated with 
each other. The ionized fraction auto-correlation func- 
tion is approximated as a Gaussian with rms given by 
the (fixed) characteristic H II region size. Under these 
assumptions the kSZ power spectrum can be calculated 
analytically. The quantitative details of the signal pre- 
dicted by this model depend on the assumptions made 
about the typical ionized patch size, and the time and 
duration of reionization, but generically it predicts a dis- 
tribution which is very strongly peaked, more so than our 
simulations find. This is related to the assumed Gaus- 
sian distribution around a fixed characteristic size of the 
ionized patches. The actual simulations also find a char- 
acteristic size for the ionized bubbles, but one that is 
also evolving with redshift, and different size distribu- 
tions around it, which yields somewhat less sharp peaks 
t han this analytical m odel predicts. 

ISantos et alJ l)2003|) proposed another simple model for 
the reionization patchiness, which assumes that the cor- 
relation function of the ionized patches is proportional to 
the one for the density field, with a bias factor which is 
time-dependent, but not scale-dependent. The power is 
filtered at small scales, below the typical size of the ion- 
ized bubbles using a Gaussian filter in k-space. The time- 
dependence of the typical size of the ionized patches is as- 
sumed to be proportional to [1 — x m (i)]~ 1 ^ 3 , where x m (t) 
is the mean mass- weighted ionized fraction at time t. The 
last function i s based on the semi-anal ytical models of 
reionization of lHaiman fc Holderl |2£)03). The resulting 
kSZ power spectra are fairly flat, with much less well- 
defined characteristic scale than our simulation results, 
as should be expected based on the quickly-evolving typ- 
ical patch size assumed in this model. The peak ampli- 
tudes of the results is somewhat higher than the ones we 
find, by factors of ~ 2 — 3, probably due to their assumed 
values of the eff ective bias. 

More recently. McQ uinn et alj_([ 2005') ap plied the semi- 
analytical reionization model of Furlanetto et al. ( 2004b) 
to estimate the kSZ signal. They found lower fluctua- 



tions in their "single reionization episode" models (by 
factor of — 3), although these models have similar mean 
reionization histories (i.e. evolution of the mean volume- 
weighted ionized fraction) to our simulations. Their 
results also predict a peak at somewhat larger scales 
(£ -2000) than ours (I -2000-9000). Their "extended 
reionization" scenarios peak to similar values to the ones 
wc find ([£(£ + l)Q/(27r)] max - 10~ 12 ) again at I -2000, 
but have a very different assumed reionization histories 
than any of our simulations. They found that the kSZ 
post-reionization contribution dominates the patchy sig- 
nal at the scales of interest, while we find that the two 
signals are comparable in magnitude. 

6. SUMMARY 

We have derived the kSZ CMB anisotropics due to 
the inhomogeneous reionization of the universe. This 
is the first such calculation based on detailed, large- 
scale radiative transfer simulations of this epoch. As 
we have shown above, these simulations follow large 
enough volume to capture the full range of scales rele- 
vant to the large-scale reionization geometry. They also 
include most of the velocity fluctuations power, which 
is not the case for smaller-box simulations. We have 
also approximately corrected for the velocity power still 
missing from the box. The resulting sky power spec- 
tra peak at £ = 2000 — 8000 with maximum values of 
[£(£+l)Ci/(2Tr)] max - lxl0~ 12 . The angular scale of the 
peak roughly corresponds to the typical ionized bubble 
sizes observed in our simulations, which is — 5 — 20 Mpc, 
depending on the assumed source efficiencies and the gas 
clumping at small scales. The kSZ anisotropy signal from 
reionization dominates the primary CMB signal above 
£ = 3000. At large scales the patchy kSZ signal depends 
largely on the source efficiencies and the large-scale ve- 
locity fields. It is higher when sources are more efficient 
at producing ionizing photons, since such sources pro- 
duce large ionized regions, on average, than less efficient 
sources. The introduction of sub-grid gas clumping in 
the radiative transfer simulations produce significantly 
more power at small scales, but has little effect at large 
scales. The sub-grid gas clumping also significantly en- 
hances the non-Gaussianity of the kSZ maps, resulting in 
up to one order of magnitude more of the brightest kSZ 
regions than a Gaussian would predict. The integrated 
kSZ signal is strong enough to be detected by upcoming 
experiments, like ACT and SPT. However, the separation 
of the patchy reionization signal from the contribution by 
the fully-ionized gas after reionization seems difficult. 

Our current simulations do not include the smallest 
atomically-cooling ionizing sources, these with total mass 
10 s M Q <; M tot < 2.5 x 10 9 M Q , as these are not resolved 
in our base N-body simulations. These smaller sources 
have important effects on the early global reionization 
history, but are also strongly suppre ssed due to Jeans - 
mass filtering in the ionized regions ijlliev et al.l l2~006c). 
During most of the evolution the reionization process is 
thus dominated by the larger sources, which are resolved 
here, and which dictate the large-scale reionization geom- 
etry. Consequently, we do not expect that the presence of 
the low-mass ionizing sources would change our current 
conclusions significantly. 
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